NECROTIC TUMOR GROWTH: AN ANALYTIC APPROACH 



MARTIN KOHLMANN 

Abstract. The present paper deals with a free boundary problem modeling 
the growth process of necrotic multi-layer tumors. We prove the existence of 
flat stationary solutions and determine the linearization of our model at such 
an equilibrium. Finally, we compute the solutions of the stationary linearized 
problem and comment on bifurcation. 
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1. Introduction 

Mathematical models for tumor growth have been considered with regularity the 
the applied sciences literature in recent years. From the mathematical point of view, 
free boundary models are of particular interest: In these models, a tumor cell at 
time t > is identified with an open domain f2(t) C K n , for some n > 1, with initial 
configuration Q(0) = CIq. For simplicity, it is assumed in many models that the 
growth process of the tumor is controlled by only two quantities: the concentration 
of nutrient (e.g., glucose or oxygen), denoted as a(x,t), and an internal pressure 
p(x,t), which both have to solve an elliptic problem on the time-dependent and 
unknown domain fi(t), with suitable conditions on the free boundary dfl(t). Finally, 
an evolution equation for the free boundary dQ(t) is needed, and usually it is 
derived from a simple application of Darcy's law pertaining to the fact that the 
tumor behaves as an incompressible ideal fluid. 

In many publications dealing with free boundary problems for tumor growth, 
the domain £l(t) is assumed to be spherically symmetric and n = 1, cf. the seminal 
papers [H, S, E3- The present work is innovative for the following three reasons: 
• We are looking at strip-shaped tumors: Lately, biologists have discovered 
that a peculiar kind of in vitro tumors can by cultivated by a special tissue 
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culture technique, cf. |12Hl4l |. In our model, we imagine the tumor to start 
growing from the flat bottom of a Petri dish. 

• An additional feature of our model is that we distinguish between a necrotic 
core, localized at the bottom of the Petri dish, and a non-necrotic shell 
which is lying above. In consequence, our problem has two free boundaries 
confining a time-dependent domain on which we study elliptic problems for 
nutrient and pressure. 

• Finally, we present a two-dimensional model (i.e., n = 2). 

We refer the reader to where the authors explain a sophisticated approach to 

the growth of non-necrotic multi- layer tumors, and where spherically symmetric 
necrotic tumor cells are studied. Based on the model assumptions in [3|, Iff, we now 
present the following problem: 

Let § = R/27rZ and consider two positive time-dependent functions p\ < p 2 on 
S. Let furthermore 

n w ,/*(*) = {(%,y) G M 2 ; x G S, pi(t,x) <y< p 2 {t 1 x)} 

with the boundary components 

r w (t) = {(x,y) e K 2 ; y = Pi {t,x)}, i = 1,2. 

The outward unit normal of T Pi (t) with respect to Q Pl iP2 (t) is denoted by i>i , for 
i = l,2. We obtain v\ and v 2 by computing the gradients of the functions Ni(x, y) = 
V- Pi{x): 

VTVi , ViV 2 
v\ = — . and v 2 



|VJVi| [VJV 2 | " 

We will write m = — VTVi and n 2 = \7N 2 . Let m denote the curvature of T Pi (t). 
It is well known that Ki can be computed explicitly using the formula 

_ _ Pixx • _ i 9 

(i + pL) 3/2 ' ' 

Let < pi.o < p 2 ,o be periodic functions on R so that pi(x,0) = pi,o(x). The 
nutrient a should satisfy a stationary diffusion equation. Furthermore, we assume 
that there is a constant supply a > of nutrient on T P2 (t) and that the normal de- 
rivative of a vanishes on T pi (t). Next, the Laplacian of the pressure is proportional 
to the difference a — a, with proportionality factor — p,; here a and p are positive 
parameters. The reason for this assumption is that if a < a, then the tumor volume 
locally decreases, whereas the tumor grows in regions where a > a. The boundary 
conditions for the pressure are the so-called Laplace- Young conditions: We assume 
that the pressure on T P2 (t) is proportional to k 2 ; the proportionality constant is 
the surface tension coefficient 72 > 0. Similarly, we have p = J1K1 — c, with 71 > 
and c a positive constant. We also assume that 

(1) a>a 

to have a reasonable long-time behavior, as explained in Q. Finally, the normal 
velocity of the boundary components is equal to the cell movement velocity in 
the direction n 2 on T P2 and — ni on T Pl respectively. This yields two evolution 
equations for the moving boundaries. 
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Figure 1: A free boundary problem for the growth of multi-layer tumors with a 
necrotic core. 



Our mathematical model is given by the following system of equations: 
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A solution to ([2]) is a tupel (ct(x, y, t),p(x, y, t), pi(x, t), p 2 {x, i)), where t € [0, T], 
T > 0, x € § and pi(x 7 t) < y < p 2 {x,t), so that all equations of ((5J) are satisfied 
pointwise. In the following sections, we will only discuss the stationary version of 
(J2J. It is given by the following system of equations 

Act = ct 

-fi(a - ct) 



c 



(3) 
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dm 
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r 

1 P 2 ! 



7i«i 

CT 


72^2 

where (a,p,p 1: p 2 ) = (a(x,y),p(x,y), px(x), p 2 (x)). 

The layout of this paper is as follows: We first prove the existence of fiat station- 
ary solutions, i.e., solutions (a,p, pi, p 2 ) with ct = cr(y), p = p(y) and with constants 
< p\ < pi- Precisely, for any given set of positive values ct > ct and // there is 
a flat stationary solution which is unique up to a shift in the y-direction. Next, 
we linearize the system ((3]) at such an equilibrium and use Fourier expansions to 
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obtain solutions of the linearized system. The calculations carried out in the fol- 
lowing sections especially show how to choose c and the surface tension coefficients 
71 and 72 to obtain non-trivial solutions. In an outlook, we present the bifurcation 
problem associated with the model ([3]). 



Acknowledgement. The author thanks the anonymous referee for asking about 
the bifurcation problem associated with the model of the paper at hand which led 
to an additional chapter compared to the initially submitted version. 



2. Flat stationary solutions 

Let (<7*,jj*, pi*, p 2 *) be a flat stationary solution of the problem i.e., we have 
that 

c* = CT * , 
p'l = 

C*(p2*) = 

<{pu) = 0, 
P*(pi*) = -c, 



P*{P2*) 
P'*(Pl*) 

I p'M 



We first solve the subproblem for the nutrient concentration and find that 

cosh y — tanh pi* sinh y 



(4) 



coshp2* — tanh pi* sinhp2* 



is the unique solution of 



Next we consider the boundary value problem 





= (7* 


(T*{P2*) 


= CT, 


<(P1*) 


= 0. 



p'l 



= -p(cr* - cr), 
P*(P2*) = 0, 

p'*{pi*) = 0, 



p*(y) = pi' 7 - CT *) + p°pi*(p2* -y) + \pv{ y 2 - pjj. 



which has the unique solution 
(5) 

Since we must demand p' se (p2*) — 0, we get the condition 

a tanh(p 2 * - pi*) 

(6) - = ■ 

<J P2* — Pi* 

Letting P2* — pi* = 6, it follows from < a < a that Eq. ^ has a unique solution 
5 € (0, 00). Next the constraint p*(pi*) = — c results in the condition 

1 _ r2 _ cosh pi* — tanh pi* sinh pi* 

(7) pa — -pad = per —. -7 : ^-—. ft — c - 

z cosh(px» + 0) — tanhp!* sinh^pi* + 0) 
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In view of the addition theorems for hyperbolic functions and the relation ([6]). 
Eq. ([7]) can be simplified to 



(8) 



pa 
cosh 8 



1 — cosh 8 



1 



- 8 sinh 8 



Since the term in brackets is positive for any 8 > 0, we have c > 0. Moreover, there 
is no condition on pi*, so that we obtain a flat stationary solution of @ for any 
fixed pi* > 0. This provides a proof of the following theorem. 

Theorem 1. Fix p, a, <r, 71, 72 > and assume that ((T|) holds true. Let 8 be the 
solution of * = ta " h s and define c according to Eq. l[8]), TTien there is a one- 
parameter family of fiat stationary solutions px*, pi* +8) to Eq. ([2]), where 
<r* and p* are given by Eq. (j4]) and Eq. ([5]) ; respectively. 

3. The linearized problem and its solutions 

A standard technique to tackle moving boundary problems is to transform the 
problem under consideration to a problem on a fixed (and preferably simple) ref- 
erence domain, to solve the problem on the reference domain and to transform its 
solutions back to obtain solutions of the original problem. 



Assume that px,p2 € C+(8) and let y) 



so that the map 



ip: (x,y) H> (x',y') := (x,"&(x,y)) establishes a C 2 -diffeomorphism O pliP2 -4§x 
(0, 1). The strip ft = § x (0, 1) will be our reference domain, with the boundary 
components Ti = § x {0} and T2 
and that ip(T Pi ) = I\, for i = 1, 2. 



p2 



^Pl,P2 



Pi 



= § X {1}. 


In particular, 


we have that I 1 , ~ 




\ V ' 








r 2 






n 




X 




\ r i 



Figure 2: Transformation of the tumor domain onto a fixed reference domain. 



Let us further introduce the operators 

A(pi,p2)u — [A(uo ip)] o ip' 1 and Bi(pi, pi)u = (tr*[V(u o ip)],rii} o ip^ 1 , 

where tr* denotes the trace with respect to T Pi , for i = 1,2, and u £ C 2 (fl). A 
straightforward computation shows that 

A(pi,p2)u = u X ' X i—2u x / y 



y'{p' 2 -p'x) + p'i , .. i + [y'(p'2-p'i) + p'i? 



pi - pi 



+ Uylyl - 



a / 2 "L h/'frWO 



(P2 - Pi) 2 



Pi] 



(P2 - Pi) 2 



P2 - Pi 
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and 



B 1 (p 1 ,p 2 )u = p'iu x \ y , =0 - 

B 2 (pl,P 2 )u = -p 2 Ux'\yi = l 

Hence the transformed problem reads 



P? 



1 



Pi - Pi 
P2 - Pi 



y> |l/'=0) 



lyi \y'=l- 



(9) 



A{pi,p 2 )<y 


= a 


in 


tlx [0,2], 


M,p\,P2)v 


= -P(° - o-) 


in 


fix [0,7], 


B 1 (p 1 ,p 2 )(T 


= 


on 


Ti x [0,7], 


p 


= 71/ci-c 


on 


Ti x [0,7], 


a 


= a 


on 


r 2 x [o,7], 


P 


= 72^2 


on 


r 2 x [o,j], 


Pi,t 


= I3 1 (p 1 ,p 2 )p 


on 


Ti x [0,7], 


P2,t 


= -B 2 (p!,p 2 )p 


on 


r 2 x [o,7], 


Pi 


= Pl,0 


for 


t = o, 


P2 


= P2S1 


for 


t = 0. 



We now pick a flat stationary solution (<7*,p*, pi*, p 2 *) &s obtained in Theorem Q] 
and, for e > 0, we let 



(10) 



/ cr(x',y',t) \ 
p(x',y',t) 

pi(x',t) 
V p 2 {x',t) J 



I a*(y'S + Pl *) \ t X(x',y',t) \ 
P*(y'$ + Pi*) +£ P(x',y',t) 
Pi* r{x',t) 

V P2 * J \ s(x',t) J 



where S = p 2 * — p\* > 0. Our regularity assumption on the new unknowns 
(£,P,r,s) is that £,P g C([0,T]; C* 2 (ft)) and r, s G C^^T]; C 2 (§)). We also 
introduce the second order linear differential operator b TtS on C 2 ([0, 1]) given by 

M«) = l( s - r)«'V* + Pi*) + ll/V' - O + r>V* + pi*). 
o 

The linearization of the problem (|9]) at the flat stationary solution (cr^p*, p\*, p 2 *) is 
obtained by inserting (|10|) into the problem ([9]) and by differentiating each equation 
with respect to e at e = 0. This yields 



(11) 



1 V , , 

Jl^y'y' 



& r , s (CT*) + £, inOx[0,T], 
Z y ,(x',0,t) = E(x',l,i) =0, onSx [0,T], 

iW + £iW = &r, s (p*) - MS, in il x [0,T], 
P(>',0,t) = -^xr"{x',t), on S x [0,T], 
P(a;',l,i) = --y 2 s"(x\t), onSx [0,T], 



r t (a:',t) 

S t (x',t) ■ 

r(x',0) = r (x'), on § 
s(x', 0) = SqI^Oj on S 



\P y .{x',0,t), onSx[0,T], 
$P y ,(x',l,t), onSx[0,T], 
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The stationary version of (fTTj) is 

^x'x' + Jz'Ey'y' = b r , s ((T*) + S, ill fi, 

Xy(x',0) = £(x',l) = 0, on S, 



(12) 



Px'x' H~ gs-Ey'v' — ^r,s(p*) M^J m ^) 
P(x',0) = -7i r", onS, 
P(x', 1) = -72s", on S, 



P yf (x',0) = 0, on§, 
P^Oe'.I) = 0, onS. 

To obtain solutions of (H2), we expand (£, P, r, s) as Fourier series and denote the 
coefficients with respect to the basis functions {1, cos(fcn'), sin(fcn'); k £ N} by 

{A (y'), A k (y'), B k {y')- k £ N} for the variable S, 
{Mo(y'), MfcCtf'), iVfc(vO; fc e N ) for thc variable P, 
{ao,afc,6fc; fc € N} for thc variable r, 

{co,c k ,d k ; k £ N} for the variable s. 

We now proceed in the following steps: First, we solve the boundary value problems 



Am 

A k (l) 




A k (y') 

0, 

0. 

B k (y>) 

0, 

0, 



9k(y'), 



for fc = 1,2, 



here 



A-(y') = ^(cfe - a k )a'l(y'8 + p x .) - k 2 [y'{c k - a k ) + a fc ]<(y'<5 + p x .) 



2cr(c fc - Ofc) 



cosh(y'<5) 



k 2 aa k 



sinh(y'<5) 



fc 2 cr(c fe - a fe ) , 



y' sinh(y'(5), 



5 cosh 5 cosh J cosh 5 

where we have used once again the addition theorems for hyperbolic functions. Thc 
9k(y') are similar; we simply have to replace the a k with the b k and the c k with thc 
d k . It is straightforward to obtain that 

ct(cq - ao) 1 



Mv') 



cosh S 



■ [y' siah(y'S) — tanh<5cosh(y'<5)] 



and 



A k {y') = 



aa k tanh((5Vl + k 2 ) ac k tanh S 



y/1 + k 2 cosh S cosh(<5vT+~P, 

— ^ vTOfeV) + -^j 

coshd 



(15) 



vT + A; 2 cosh <5 

cosh * y smHyS) - 



cosh(<5\/l + fcV) 
sinh(j/'<5) 
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The Bk emerge from the Ak by exchanging the a k with the b k and the c k with the 
d k . We next turn our attention to Mq which has to satisfy the following conditions: 



(16) 







A*o(0) 


= o, 


M^(l) 


= o, 


M (0) 


= o, 


, M (l) 


= 0. 



cosh 5 



-)• 



Explicit calculations show that there is a solution Mq if and only if cq = a,Q and 
precisely Mq = 0. Next, we plug the solutions A k and B k into the problems 



(17) 
and 
(18) 

for k = 1,2,. 
/fc(y') = 



-fc 2 Af fc (y') + £M£(|/0 

K(0) 
K(l) 



-fc 2 W) + ^7V£(y') 
^(1) 



o. 

0. 



o, 

0. 



<?fc(z/'), 



here 



2 1 1 



(c k - a k ) (a- - a*(y'5 + p u )) 

^k 2 (y'(c k - Ofc) + a k ) {aSy' - a'^y'8 + p u )) 



fj,k (TOk . 



-zpa(c k - o fc ) - pk 2 aSa k y' - pk 2 aS(c k - a fc )i/ 2 + 
d cosh o 

2pa{c k - a fe ) cosh(y , 5) + MfcMcfc-afc) ^ s . nh(y/(5) 



sinh(y'<5) 



(5 cosh <5 cosh £ 

and g k (y') is obtained as before. Again, it is straightforward to derive the solutions 



AW) 



/i cosh(y'Jfc) 
<5fc sinh(<5fc) 



er<5a fc 



cosh (5 cosh((5Vl + fc 2 ) 



cr&Xfc cosh((5fc) 



+ercfc<5'\/ 1 + k 2 tanh(5tanh(<5-\/l + fc 2 ) + crcfe tanh<5 — <5crcfc 



/icafc 



■ sm 



VI + fc 2 cosh 5 

Ck tanh <5 



h(y'Sy/l + k 2 ) 



ak tanh((5v 1 + k 2 
3 h(oVl + P) Vl + fc 2 coshT 



3h( 2 /'5vT 



(19) 



sinhfy ok) smh(oy ) 

k cosh 5 

+p,aSa k y' + paS(c k - a k )y' 2 , 



pa(c k - a k ) , . , 

coshj y smh{Sy } 



and N k accordingly. Since we are looking for non-trivial solutions, we assume that 
(afc,&fc), (cfc,dfc) 7^ (0,0) and that akdk = bkCk- Without loss of generality, we will 
henceforth suppose that a k ^ and Ck ^ 0. Next we may choose the parameters 
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71 and 72 so that 



7i 



k 2 



« A: 



-M fc (0) and 72 



A' 2 



-M fc (l); 



precisely, 

7i : 



c>fc 3 sinh((5fc) 



+ 1 



cosh((5fc) tanh<5 



Sk tanh((5Vl + fc 2 ) sinh(Jfc) 
Vl + fc 2 cosh (5 



cr , fcsinh((5fc) -VT+ fc 2 sinh(cVl + fc 2 ) 

— + tanhd ^=^= 

cr cosh(<5v 1 + fc 2 ) 



(20) 
and 



72 



fia 



cosh d cos h^A/ITX 2 ) 
afc / tanh<5 



Sk 3 tanh(<5fc) 



(21) 



c fe \cosh(5k) cosh5cosh(5Vl + k 2 ) 
+<5 tanhd 7 (fctanh^fc) - \fl + k 2 tanh(cVl + fc 2 )) + 5 - tanhtf 



For simplicity, we will now furthermore assume that = for any k £ N. It 
is easy to see that the terms in brackets standing one underneath the other in 
the formula for 71 converge to +00, S — (1 + 4^-) tanh<5 and zero, respectively, for 
fc — > 00. Similarly, one deduces that the terms in brackets for 72 tend to 5 — tanh5 
as k approaches 00. It follows that 71,72 > for k sufficient large. Moreover we 
have that 71, 72 — > 0, as k — > 00. The figure below shows the functions 71 (fc) and 
72 (fc) for a particular choice of the parameters a, a and /x. 




10 



12 14 



Figure 3: The coefficients 71 and 72 (dashed) for a = g<7 = fx = 1. 

Note that we have J^ fe |afc| 2 < 00 and \ck\ 2 < 00, so that afc and Ck are null 
sequences, but in general, it is difficult to obtain control of the ratios ^ and 24, 
which is why we are working with the additional assumption = afc here. 

Our result can be formulated as follows. 
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Theorem 2. Pick fceN and let ak,bk,Ck,dk be real constants satisfying the rela- 
tions 

ai + bl^0, cl + dl^0, a k dk = bkc k . 
The linearized problem (|12p has a nontrivial solution (£, P, r, s) if and only if 71 
and 72 are as in (|20p and (|2ip and i/iis solution is given by 

E(x', J/') = Afe(y') cos(fca;') + sin(fcx'), 
F(x' ;2 /') = M k {y')cos(kx') + N k (y') sm(kx'), 

r(x') = afc cos(fcx') + b k sm(kx'), 

s(x') = c k cos(kx') + d k sin(fcx'); 

the coefficients A k and M k are as in (| 1 5|) and (|19p and and iV^ are obtained 
by exchanging a k with b k and c k with d k . If a k = Ck ^= or b k = dk 7^ and k is 
sufficiently large, we have 71,72 > and 71,72 — > 0, /or fc — > 00. 

4. Outlook: The bifurcation problem 
The results of Theorem [2] motivate to study bifurcation for the problem ([2]) ; 



cf. [16[ where bifurcation for the non-necrotic version of our strip-shaped tumor 
growth model is established. In some older papers, bifurcation solutions for radially 
symmetric free boundary value problems have been constructed by using a power 
series technique; see, e.g., The modern method of analysis in (3] is based on 

an application of the following theorem of Crandall and Rabinowitz. 

Theorem 3 (see 0). Let X and Y be real Banach spaces and let G(u, A) be a C q 
map (q > 3) from a neighborhood of a point (uq, Xq) £lxl into Y . We assume 
that 

(1) G(u o ,Xo) = G x (u o ,Xo) = 0, 

(2) Ker G u {uq, Xq) is one- dimensional and spanned by uq, 

(3) Im G u (uo,Xq) has codimension 1 and 

(4) Gaa(uo,A ) e Im G u (u Q ,X ), G u \(u Q ,X Q )u £ Im G u (u 0l X ). 

Then (uq 7 Xq) is a bifurcation point of the equation G(u,X) = in the sense that 
in a neighborhood of (uo,Xq), the set of solutions of G(u, A) = consists of two 
C q ~ 2 smooth curves Ti^ which intersect only at the point (uo,Ao) and can be 
parameterized as follows: 

Ti : (u(A),A), I A — Aq | is small, u(Xo) — Uo, u'(Xq)=0, 
T 2 : (u(e), X(e)), |e| is small, (u(0),A(0)) = (u ,X ), u'(0) = u . 



In [16|, the authors let the surface tension coefficient 7 (which corresponds to 
72 in our model) play the role of the bifurcation parameter A. Our crucial problem 
is that there are two surface tension coefficients 71 7^ 72 in the necrotic variant of 
the multi-layer tumor growth model which suggests that the Crandall-Rabinowitz 



Theorem is not suitable for our purposes. Moreover, the technique presented in [16 1 



is fairly standard to obtain bifurcation branches for related free boundary models 

where 



11 



and has already been applied in various publications 0, Ej]; see also 
bifurcation from radially symmetric solutions of a necrotic tumor growth model is 
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discussed. Because of that, we provide some supplementary material concerning 
the model ([2} in this outlook and prepare a functional analytic formulation which 
might be suitable to derive bifurcation. It remains an open problem to establish 
the existence of bifurcation branches for which we probably need some deep and 
new ideas. 

Let h m+a (E>), m <G {2, 3, . . .} and a <E (0, 1), denote the little Holder space on the 
circle, i.e., the closure of C°°(S) in the Holder space C m+a (E>). The cone of positive 
functions in h m+a (S) is denoted as /i™+"(§). We define h m+a (Ti) analogously. 
The small Holder spaces are used frequently since they are Banach algebras (under 
pointwisc multiplication) and the embedding h r (S) h s (B), r > s, is compact. To 
keep our notation as simple as possible, we will label the coordinates in ft by x and 
y in the sequel. 

First, we establish that the linear second-order differential operator 

A( Pl ,p 2 ) = a n dl + 2a 12 d x d y + a 22 d 2 y + bd y , Pl ,p 2 e ^ i+Q (§), 

introduced in Section |3] is uniformly elliptic in ft as defined in [lo| . p. 30. Pick 
£ = (£1,^2) € M 2 \{0}. Since the coefficients a,ij = a,ji, for i,j = 1,2, are continuous 
functions on S x (0, 1), it is clear that . ctij£i£j < A|£| 2 , for some A > 0. On the 
other hand, we have 

f = fa - yAEk^l±AzX + ^-L^e 2 2 = Hi 2 , 

^ V P2 - Pi J (P2- Pi) 2 



where 



/ 1 y{p' 2 -Pi)+Pi \ 



\ P2-P1 / 

Clearly, A is invertible and there is A > such that 

£«^>^4^I£I 2 >A|£I 2 ; 

here, lA -1 ^ is the Frobenius norm of A -1 . 

For any given pi, p 2 € h™ +a (§) we solve the elliptic boundary value problem 

!A{pi,p 2 )cr = a in ft, 
Bi(pi,p 2 )a = on§, 
c\y=i = a on S, 

and obtain a unique solution a <G /i m+Q (ft). We let TZ denote the solution operator 
for the nutrient concentration and write a = lZ(fti, p 2 )<5- It follows from elliptic 
regularity theory [lfj that 

K(; -)a £ C 00 ^™^) 2 , h m+a (H)). 

Second, we study the Neumann problem 

A(pi,p 2 )p = -p,{TZ{pi,p 2 )a - a) in ft, 
I B 1 (p 1 ,p 2 )p = on§, 
B 2 (pi,p 2 )p = on S, 
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which is solvable if and only if 

$(Pi,P2):= / (K(pi,p 2 )o- - a) (p 2 - Pi)dxdy = 0. 
J n 

Here, $ G C°°(/i" l+Q (§) 2 , R). Using gj) and ©, we compute 

$(pu,P2*) = 5 I (H(pi*,P2*)v ~ cr) dxdy 

cosh(y<5 + pi*) — tanh pi* sinh(j/5 + pi*) tanh<5 



= 5a 



cosh^* — tanh pi* smhp 2 * 



dx dy 



= 



(24) 

and, in view of 



*'(pi.,pa.)(0,l) 



lim - [$(pi*,p 2 * +e) - $(pi*,p 2 *)] 

e-s-0 £ 



lim - / [7^(pi»,p 2 * + e)er — a] (5 + e) dxdy 
£ ,/n 

(5 /* 

lim-y [7£(pi», ,02* + e)a — a] dxdy 

I [7£(pi*,p 2 * + e)o - TZ{pi*, P2*)v] dxdy - 
I lim - [7£(pi*,p 2 * + e)a - K(pi*, p-2*)v] dxdy 



lim 



;*(pi*,p 2 *) 



all 



d 

n ds 
1 



e=0 



cosh(y(5 + e) + pi*) - tanhpi* sinh(y(<5 + e) + pi*) 
cosh(p 2 * + e) — tanhpi* sinh(p 2 * + e) 



tanh 5 — tanh" 5 



dx dy 



(25) 

the explicit calculations in the last steps are left to the reader. The expression 
obtained in Eq. (f25|) is nonzero, since 1 — i tanh 6 — tanh 2 6 = implies that 
5 = sinh<5cosh5 which is possible only if 5 = 0. 

It follows from ([25]) and the continuity of $' that there is a neighborhood U 
of (pi*, pa*) in /i™ +Q (§) 2 such that Vfa, p 2 )(0, 1) 7^ for all (pi,p 2 ) G 17. Thus 
^(Pi; P2) = defines a smooth Banach submanifold M. of codimension 1 in a small 
neighborhood of (pi*,p 2 *), i.e., 

X = |(p 1 ,p 2 )e^ +Q (S) 2 ; maxllpi-pi*!^^ < S, $(p 1 ,p a ) = o|. 

For any (pi,p 2 ) G the problem (|23p has a solution which is unique up to a 
constant. Let T(pi,P2) be the solution operator for ((23)) which associates to the 
right-hand side — / the solution p which is zero at the origin. We then have 



V = pT(pi, p 2 ) {K(pi,p 2 )a - cr) 
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and 

T e C™{h™ +a {§) 2 X{h m+a {Ti),h m+2+a {Ti)). 

Let trj, i = 1,2, denote the trace operator for Ti and r 2 respectively. We now set 

Si{pi,P2) = fJbtXi o T(/9i,/3 2 ) (lZ(p ll p 2 )d- - a) 
and recall that the curvature operator is given by 



K(p) 



cPp 
dx 2 



dp 
dx 



-3/2 



We have shown that the problem © can be rewritten as 
(26) 



S\{pi,P2) 
S2{Pl,P2) 



Cq + C 
CO 



7i«(Pi), 

l2K(p 2 ), 

0, 



where £ C°°{M, h m+a (S)) and k e C°°(7W, /i m - 2+a (§)). With Si 
and S-2 = S2 + cq, we conclude that (|26p is equivalent to 



Si 



- cq 



(27) 



F{pi,p 2 , 71=72) 



5i(/?i,p 2 ) - 7i«(pi) 

^2(Pl,P2) - 72«(P2) 



and F e C°°(7W x E^/i 15 



'(S) 2 ) 



16j can be rewritten as the zero level set of a 



While the non-necrotic model in 
function X x R + —> Y, where X and Y are real Banach spaces suitable for the ap- 
plication of the Crandall-Rabinowitz Theorem, the mapping F in (|27[) employs two 
positive parameters 71 , 72 as inputs which is not compatible with the assumptions 
of Crandall-Rabinowitz. Thus the bifurcation problem for (|2T|) remains a subject 
for further research. 



References 

[1] Byrne, H., Chaplain, M.: Growth of nonnecrotic tumors in the presence 
and absence of inhibitors. Math. Biosci. 130, 151-181 (1995) 

[2] Crandall, M.G., Rabinowitz, P.H.: Bifurcation from simple eigenval- 
ues. J. Funct. Anal. 8 321-340 (1971) 

[3] Cui, S., Escher, J.: Wcll-posedncss and stability of a multi-dimensional 
tumor growth model. Arch. Rational Mech. Anal. 191, 173-193 (2009) 

[4] Escher, J., Matioc, A.: Bifurcation analysis for a free boundary prob- 
lem modeling tumor growth. Arch. Math. 97 79-90 (2011) 

[5] Escher, J., Matioc, A., Matioc, B.: Analysis of a mathematical model 
describing necrotic tumor growth. arXiv:1005.2506vl [math.AP] 

[6] Friedman, A., Hu, B.: Bifurcation for a free boundary problem mod- 
eling tumor growth by Stokes equation. SIAM J. Math. Anal. 39(1) 
174-194 (2007) 

[7] Friedman, A., Hu, B., Velazquez, J.J.L.: A Stefan problem for a pro- 
tocell model with symmetry-breaking bifurcation of analytic solutions. 
Interfaces Free Bound. 3, 143-199 (2001) 



Necrotic tumor growth: an analytic approach 



14 



[8] Friedman, A., Reitich, F.: Analysis of a mathematical model for the 
growth of tumors. J. Math. Biol. 38, 262-284 (1999) 

[9] Friedman, A., Reitich, F.: Symmetry-breaking bifurcation of analytic 
solutions to free boundary problems. Trans. Amer. Math. Soc. 353, 
1587-1634 (2000) 

[10] Gilbarg, D., Trudinger, N.S.: Elliptic Partial Differential Equations of 
Second Order. Springer, New York, 1977 

[11] Hao, W., Hauenstein, J., Hu, B., Liu, Y., Sommcse, A., Zhang, Y.: 
Bifurcation for a free boundary problem modeling the growth of a 
tumor with a necrotic core. Submitted to Nonlinear Analysis: Real 
World Applications (January 19, 2011). 

[12] Kim, J.B., Stein, R., O'Harc, M.J.: Three-dimensional in vitro tis- 
sue culture models for breast cancer — a review. Breast Cancer Res. 
Treat. 149, 1-11 (2004) 

[13] Kyle, A.H., Chan, C.T.O., Minchinton, A.I.: Characterization of three- 
dimensional tissue cultures using electrical impedance spectroscopy. 
Biophys. J. 76, 2640-2648 (1999) 

[14] Miillcr-Klicscr, W.: Three-dimensional cell cultures: from molecular 
mechanisms to clinical applications. Am. J. Cell Physiol. 273, 1109- 
1123 (1997) 

[15] Ward, J., King, J.: Mathematical modelling of avascular-tumour 
growth. IMA J. Math. Appl. Med. Biol. 14, 39-69 (1997) 

[16] Zhou, F., Eschcr, J., Cui, S.: Bifurcation for a free boundary problem 
with surface tension modeling the growth of multi-layer tumors. J. 
Math. Anal. Appl. 337, 443-457 (2008) 



Peter L. Reichertz Institute for Medical Informatics, University of Braunschweig, 
D-38106 Braunschweig, Germany 

E-mail address: martin.kohlmannaplri.de 



